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^^ , We present in detail a new systematic method which can be used to automatically eliminate 

00 ' the renormalization scheme and scale ambiguities in perturbative QCD predictions at all orders. 

^Nj ' We show that all of the nonconformal /3-dependent terms in a QCD perturbative series can be 

readily identified by generalizing the conventional renormalization schemes based on dimensional 

regularization. We then demonstrate that the nonconformal series of pQCD at any order can be 

resummed systematically into the scale of the QCD coupling in a unique and unambiguous way due 

to a special degeneracy of the /3 terms in the series. The resummation follows from the principal 

Oj, of maximum conformality (PMC) and assigns a unique scale for the running coupling at each 

D ' perturbative order. The final result is independent of the initial choices of renormalization scheme 

rS^ ' and scale, in accordance with the principles of the renormalization group, and thus eliminates an 

unnecessary source of systematic error in physical predictions. We exhibit several examples known to 

order a*; i.e. i) the electron-positron annihilation into hadrons, ii) the tau-lepton decay to hadrons, 

iii) the Bjorken and Gross-Llewellyn Smith (GLS) sum rules, and iv) the static quark potential. We 

show that the final series of the first three cases are all given in terms of the anomalous dimension 

p^ . of the gluon field, in accordance with conformality, and with all non-conformal properties encoded 

vQ ' in the running coupling. The final expressions for the Bjorken and GLS sum rules directly lead to 

^^ I the generalized Crewther relations, exposing another relevant feature of conformality. The static 

• . quark potential shows that PMC scale setting in the Abelian limit is to all orders consistent with 

'b! ' QED scale setting. Finally, we demonstrate that the method applies to any renormalization scheme 

and can be used to derive commensurate scale relations between measurable effective charges, which 

provide non-trivial tests of QCD to high precision. This work extends ELM scale setting to any 

perturbative order, with no ambiguities in identifying /?-terms in pQCD, demonstrating that ELM 

scale setting follows from a principle of maximum conformality. 
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I. INTRODUCTION 

An important goal in high energy physics is to make 
perturbative QCD (pQCD) predictions as precise as pos- 
sible, not only to test QCD itself, but also to expose new 
physics beyond the standard model. Recently, we showed 
a systematic method to determine the argument of the 
running coupling order by order in pQCD, and in a way 
that can be readily automatized [1, 2]. The new method 
satisfies all of the principles of the renormalization group 
[3] , and it eliminates an unnecessary source of systematic 
error. The resulting predictions for physical processes are 
independent of theoretical conventions such as the choice 
of renormalization scheme and the initial choice of renor- 
malization scale. The resulting scales also determine the 
effective number of quark flavors at each order of pertur- 
bation theory. The method can be applied to processes 
with multiple physical scales and is consistent with QED 
scale setting. 

In this paper we review the method in detail and pro- 
vide the complete generalization. We show several ex- 
amples, based on observables recently published in the 
literature to four-loop order in perturbation theory. Fi- 
nally, we demonstrate that the method applies to any 
renormalization scheme and can be used to derive com- 
mensurate scale relations between measurable effective 
charges [4-7] . The method extends the Brodsky-Lepage- 
Mackenzie (BLM) method [8] to any perturbative order 
by following the Principle of Maximum Conformality [9- 
11], without leaving any ambiguity in identifying /3-terms 
at any order in pQCD. 

Previous attempts made in the literature to extend 
BLM scale setting to higher orders [12-22], have mostly 
focused on improving convergence of the perturbative se- 
ries and not removing renormalization scheme and scale 
ambiguities. Therefore they do not in general satisfy 
the self-consistency requirements of the renormalization 
group [3], nor the initial renormalization scale depen- 
dence, which must be the prerequisite of any scale-setting 
method. The main contribution of this work is to provide 
the systematic method to eliminate the renormalization 
scale and scheme ambiguities to all orders in pQCD. 

Other recent proposals similar in spirit are suggested 
inRefs. [23, 24]. 

We start our analysis in Sec. II by introducing a 
generalization of the conventional schemes used in di- 
mensional regularization in which a constant ^6 is 
subtracted in addition to the standard subtraction 
ln47r — 7£ of the MS-scheme. The (S-subtraction defines 
an infinite set of renormalization schemes which we call 
(5-7?.enormalization (TZs) schemes; since physical results 
cannot depend on the choice of scheme, predictions must 
be independent of S. As will be described in Sec. Ill, 
the 7?,5-scheme exposes the general pattern of noncon- 
formal {/3i}-terms, and it reveals a special degeneracy of 
the terms in the perturbative coefficients which allows us 
to resum the perturbative series. The resummed series 
matches the conformal series, which is itself free of any 



scheme and scale ambiguities as well as being free of a 
divergent renormalon series. It is the final expression one 
should use for physical predictions. It also makes it pos- 
sible to setup an algorithm for automatically computing 
the conformal series and setting the effective scales for 
the coupling at each perturbative order. 

In Sec. IV we provide several examples, based on ob- 
servables recently published in the literature to order a^; 
i.e. i) the electron-positron annihilation into hadrons, ii) 
the tau-lepton decay to hadrons, iii) the Bjorken and 
Gross-Llewellyn Smith (GLS) sum rules, and iv) the 
static quark potential. We show explicitly that the fi- 
nal series of the first three cases are all given in terms of 
the anomalous dimension of the gluon field, in accordance 
with conformality, and with all non-conformal properties 
encoded in the running coupling. Moreover, the final ex- 
pressions for the Bjorken and GLS sum rules directly lead 
to the generalized Crewther relations, exposing another 
relevant feature of conformality. The static quark poten- 
tial furthermore provides an example of how the method 
can be automatized to give the PMC prediction directly 
from the number of quark flavor dependence of the ini- 
tial expression. From this example we also demonstrate 
that the PMC prediction in the Abelian limit is consis- 
tent with QED scale setting. Finally, we demonstrate in 
Sec. V that the method applies to any renormalization 
scheme and can be used to derive commensurate scale re- 
lations between measurable effective charges, which pro- 
vide non-trivial tests of QCD to high precision. 



II. THE ^-7^ENORMALIZATION SCHEME 

In dimensional regularization logarithmically divergent 
integrals are regularized by computing them in d ~ 4 — 2e 
dimensions [25-28]. This requires the following transfor- 
mation of the integration measure and introduction of an 
arbitrary mass scale fi: 



dp -^ n 



i4-2£ 



(1) 



Divergences are then separated as 1/e poles and can be 
absorbed into redeflnitions of the couplings. The choice 
of subtraction procedure is known as the renormalization 
scheme and is chosen at the theorist's convenience. To 
avoid dealing with coupling constants changing dimen- 
sionality as a function of e one rescales the couplings as 
well with the mass scale /i in the d = 4 — 2e theory. In 
particular, for QCD one rewrites the bare gauge coupling 
oo = ao/AiT = g'^ /{AiiY as: 



ao 



2e rr 



(2) 



where as is the renormalized gauge coupling under a spe- 
cific renormalization scheme S and Za^ is the renormal- 
ization constant of the coupling. The mass scale /i is 
now understood as the renormalization scale. The bare 



coupling must be independent of the arbitrary scale /i, 
thus 



^dao 
Using this and the expansions 
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it is easily derived that: 
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and the /3i coefficients are known up to /^a, or four loops 
[29]. The coefficients Pi are renormalization-scheme de- 
pendent; however, it is easy to demonstrate by a general 
scheme-transformation that the first two coefficients Pq 
and Pi are universal for all mass-independent renormal- 
ization schemes. 

In the minimal subtraction (MS) scheme [30] one ab- 
sorbs the 1/e poles appearing in loop integrals which 
come in powers of 
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where c is the finite part of the integral. Since anything 
can be hidden into infinity, one can subtract any finite 
part as well with the pole. This is equivalent to redefining 
the arbitrary scale /i in Eq. (1). The MS-scheme [31] 



differs from the MS-scheme by an additional absorption 
of the term ln(47r) — 7^;, which corresponds to redefining 
/i to: 



^^ "" ^lls exp(ln47r - 7^) 



(9) 



We will generalize this by defining the 
S-TZenormalization scheme, TZs, where one absorbs 

ln(47r) — 7b — 5, i.e. 



M 



/i^ exp(ln47r - 7b - (S) 



(10) 



where 6 is an arbitrary finite number, and by appropriate 
choice will connect all MS-type schemes. In particular^: 



7^o=MS 

'^In 47r— O'K 
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The scheme-transformation between different TZs cor- 
responds simply to a displacement in their corresponding 
scales, i.e. 



In particular: 



M52 = ^i exp((S2 - Si) 



^5 = f4m °^P(^) 



(13) 



(14) 



Since all T^-^-'s are connected by scale-displacements, 
the /?-functions of a-jz^ defined in Eq. (4) are the same 
for all TZs to any order. The index S on a-jz^ is thus 
redundant and we denote it instead as a-ji. Where it will 
not be ambiguous, we will simply use a = a-jz. 

We can find a power series solution in l/ln(/i/A) for 
a by solving the renormalization group equation pertur- 
batively. It is simplest to use the extended renormaliza- 
tion group prescription [32, 33] where one works with a 
rescaled coupling and a rescaled logarithm; respectively 



Po' 



i, = ^ln(M./A) 
Pi 



The solution up to 0{LJ^) reads [2, 9]: 



d{fis) 



h^¥ 


-lnL5) + -3 [C2+C + C2 


~{2C-\yiLs + \) 


1-C3 
2 


3C2 + 5C + 3c2 - 2 - f 3C - 


5\ 
-hiLs + - In Ls 
2/ 



InLf 



^5 



O 



c [c 



-C + 3c2 - 2 



(15) 



^ Note that we have chosen MS as the reference scheme for 'R.q. 
This is done since most results today are known in this scheme; 
however there is nothing special about MS, and TZo can be rede- 
fined to be any other MS-like scheme 



where Cj = PiPl^"^ / P\ are the rescaled /3-function coeffi- 
cients and C is an arbitrary integration constant which 
in TZs is set to C = In/^p//?! in order to reproduce the 
standard Aj^ scale [2, 9, 33] . Note that we take the 
asymptotic scale A = A^^ to be the same for any TZs- 



Alternatively, one can take the scale ^ to be the same for 
any TZg, having instead different asymptotic scales Kg. 

This solution for a{iis) holds for massless QCD or as 
long as the active quark masses are below A. In the 
case where there are active quarks with masses higher 
than A, one must take the quark threshold effects into 
account when running the coupling from A to jig- This 
can be done by e.g. solving the renormalization group 
equation with an analytic /3-function that takes quark 
masses into account [34] or by using matching equations 
at each quark mass threshold when running the coupling 
to ^s [35]. 



III. OBSERVABLES IN Us 

Consider an observable in pQCD in some scheme which 
we put as the reference scheme TZq with the following 
expansion: 



Po(Q') = a(Mo)" Vr,+i(QVA^o)«(Mo)' 



fc=0 



(16) 



where /io stands for some initial renormalization scale 
and Q is the kinematic scale of the process. The full 
pQCD series is formally independent of the choice of 
the initial renormalization scale /ip, if it were possible to 
sum the entire series. However, this goal is not feasible 

I 



in practice, especially because of the k\fi^a^ renormalon 
growth of the nonconformal series. When a perturbative 
expansion is truncated at any finite order, it generally be- 
comes renormalization-scale and scheme dependent; i.e., 
dependent on theoretical conventions. This can be ex- 
posed by using the 7?.5-scheme. Since results in any 7?.^- 
are related by scale displacements, we can derive the gen- 
eral expression for p in TZs by using the displacement 
relation between couplings in any 7?,5-scheme: 

a(^o) = a(M.) + £ l^^jj^ I,.,, {-ST , (17) 

where we used ln/ig//i^ = ~5. It is useful to derive the 
general displacement relation for a(yUo)'^ for any k as an 
expansion in a up to order a^^^: 
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Inserting this expression into Eq. (16) at each order 
a{iiQ)^ we find the expression for p for an arbitrary 5 
to order a"', that is in any T^-^-scheme, to be: 



P5{Q') =riai{psT + [^2 + npori5i]a2{psT+^ + 



r-i + nPiTiSi + (71 + l)/3or2(52 + 



n{n + l){n + 2) ' 
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where pj — Q e , the initial scale is for simplicity set 



Q , and we have defined in Eq. (16) ri{l) 



to pI ^ ^ 

An artificial index was introduced on each a and d to 
keep track of which coupling each S term is associated 
with. They are not an indication of different variables; 
i.e. oi =02 = ... and 61 = 62 — ■ ■ ■■ The use of the 
artificial indices will be made clear in a moment. 

The above expression shows the scheme depen- 
dence explicitly; e.g. if TZq = MS, then choosing 
S = 6i = In Att ~ jE will give the result in the MS-scheme. 
The initial scale choice is arbitrary and is not the final 
argument of the running coupling; the final scales will be 
independent of the initial renormalization scale. 

In a conformal (or scale-invariant) theory, where 
{ft} = {0}j the S dependence vanishes in Eq. (19). 
Therefore by absorbing all {Pi} dependence into the run- 



ning coupling at each order, we obtain a final result in- 
dependent of the initial choice of scale and scheme. The 
coefficients in the final expression will thus be equal to 
those of the conformal theory. The use of TZg allows us 
to put this on rigorous grounds. From the explicit ex- 
pression in Eq. (19) it is easy to confirm that 
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OS 



-/?(« 



dps 
da 



(20) 



The scheme- invariance of the physical prediction requires 
that dps/d5 = 0. Therefore the scales in the running 
coupling must be shifted and set such that the confor- 
mal terms associated with the /?-function are removed; 
the remaining conformal terms are by definition renor- 
malization scheme independent. The numerical value for 
the prediction at finite order is then scheme independent 
as required by the renormalization group. The scheme- 



invariance criterion is a theoretical requirement of the 
renormalization group; it must be satisfied at any trun- 
cated order of the pertubative series, and is different from 
the formal statement that the all-orders expression for a 
physical observable is renormalization scale and scheme 
invariant; i.e. dp/d/.to = 0. The final series obtained cor- 
responds to the theory for which /3(a) = 0; i.e. the confor- 
mal series. This demonstrates to any order the concept of 
the principal of maximum conformality (PMC) [9], which 



J 



states that all non-conformal terms in the pcrturbative 
series must be resummed into the running coupling. 

The expression in Eq. (19) exposes the pattern of {/3i}- 
terms in the coefficients at each order. It is possible to 
infer more from Eq. (19); since there is nothing special 
about a particular value of 6, wc conclude that some of 
the coefficients of the {/3i}-tcrms arc degenerate; e.g. the 
coefficient of /3oa(Q)^ and (3ia{Q)^ can be set equal. Thus 
for any scheme, the expression for p can be put to the 
form: 



P((9^) =ri,oa{QT + h,o + nl3or2sHQ) 



n+l 



'"3,0 + nl3ir2^i + {n + l)/3o'"3, 
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where the r^^o a-re the conformal parts of the pcrturbative coefficients; i.e. r^ = r^.o + C'({/3i}). 
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The 7?.5-scheme not only illuminates the {/3i}-pattern, 
but it also exposes a special degeneracy of coefficients at 
different orders. The degenerate coefficients can them- 
selves be functions of {/3i}j hence Eq. (21) is not to be 
understood as an expansion in {/3i}, but at pattern of 
{/3i} with degenerate coefficients that must be matched. 
The artificial indices in the expansion in Eq. (19) re- 
veals how the {/3i}-tcrms must be absorbed into the run- 
ning coupling: The different J^'s keep track of the power 
of the 1/e divergence of the associated diagram at each 
loop order in the following way; the (5^a™-term indicates 



(21) 



the term associated to a diagram with i/£i+»»-"-'s di_ 
vergencc for any power p of S. Grouping the different 
(5fe-terms, one recovers in the iVc — >■ Abelian limit [36] 
the dressed skeleton expansion. Resumming the series ac- 
cording to this expansion thus correctly reproduces the 
QED limit of the observable and matches the conformal 
series with the running coupling evaluated at effective 
scales at each order. Using this information we can re- 
arrange the expression in Eq. (21) in the skeleton-like 
expansion: 



p(O') =a{Qr [n,o + n (PodiQ) + PMQf -f p2a{Qf) r^.x + -((« + \)^la{Qf + (3 + 27i)/3o/3ia(g)') r^a 
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A. Systematic All-Orders PMC scale setting 



It is easy to see from Eq. (22) that wc can rcsum all 
Tj^i terms, which come with a linear factor of /3j, to all 
orders by defining new scales Qi at each order as follows 



I 

(for simplicity, we treat the higher-power /3j terms later) : 

r-i,oa(Qi)" = ri,oa(Q)" - na(Q)"-i/3(a)r24 , 
^•2,00(^2)"+' = r2,oa(g)"+' - (n + l)a(Q)"73(a)r3,i , 
r3,oa(g3)"+' - r3,oa(Q)"+2 - (n + 2)a(Q)"+ V(a)^4,i , 

rkMQk)'' = ruMQf - fc a{Qf-^ P{a)rk+ix- (23) 



From the scale displaeement equation (17) for a it is 
straightforward to see that: 



a{Qkf = a{Q)'' + ka{Qf-^l3{a) In %+ 
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It follows from identifying Eq. (23) with (24) that to ab- 
sorb all linear Pj terms, the scales Q^. must satisfy: 
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where r^.o are the conformal coefficient and r^+i^i are 
the degenerate coefficients of linear /3j -terms. This leads 
to the self-consistency equation for Qk'- 




To leading logarithmic order (LLO) we have: 

LLO ^fc+1,1 
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(27) 



This resums all linear /3j terms, but introduces higher- 
power /3j terms beyond the order 0*^+^. For example, 
suppose that we are computing an observable to order 
aP. The scales Qk must resum all j3jr k+1.1 terms without 
introducing higher order ones up to order a^ . This means 
that Qk must be computed to the p— {k + \) logarith- 
mic order (Np~(''+^^LLO). Let us explicitly perform the 
resummation up to a^ for the ffist scale Qi, that is, up 
to next-to-next-to- leading logarithmic order (NNLLO). 
The general expression for the NLLO scale reads: 
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To find the NNLLO scale, we first write the self- 
consistency equation (exposing one higher logarithmic 
order in the denominator): 
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Then we replace the logarithms in the denominator with the expansion of its NLLO expression in Eq. (28): 
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We thus get: 
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This procedure iterates to any desired order. 
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For observables, where the higher-power /3j coefficients 
vanish, i.e. r3^2 = ?'4,2 = ^4,3 = (this is the case in e.g. 
the Adler £>- function) , these scales give the final PMC 
expression for the observable, which is invariant under 
any scheme transformation: 



P{Q^) ^ri,oaiQn,NNLLoT + r2S)a{Qn+l,N LLo) 
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This is the conformal series with coefficients that are in- 
dependent of the renormalization scheme. Note that the 



last scale remains ambiguous. This ambiguity only af- 
fects the highest order term. The final expression and 
coefficients are therefore not affected by the ambiguity 
of the last scale and thus the renormalization scheme de- 
pendence has been eliminated and the renormalization 
scale dependence only resides in the highest power cou- 
pling of the perturbative series. We note that one does 
not need the full expression of the a^ coefficient to set 
the last scale, Qn-i-s, but only the coefficient rs^i. 

Let us now generalize to observables that do depend 
on higher powers in j3j. This is for example the case in 
-Re+e--^hadrons(s)- One Can usc the procedure just de- 
scribe, but instead of Eq. (23) we use its generalization: 



rkMQk) = rk,oa{Qy - k a{QY ^/?(a)rfc+i,i + - 
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(33) 



It is easy to verify that these expressions, which define the PMC scales Qk, correctly resum all {/3i}-terms in p. Eq. (33) 
is systematically derived by replacing the In-* Q\IQ^ by r^j- in the logarithmic expansion of a{QkY in Eq. (24) up to 
the highest known rfc,„-coefiicient in pQCD. We introduce a short-hand notation of Eq. (33): 



liQk)" = a{QY + k a(Q)'=-i/3(a) {i?,a + A[}\a)Rk^2 + Af (a)i?fc,3 + • • • + A["\a)Rk,n+i} , 
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(35a) 
(35b) 

(35c) 



Following the same procedure as before, one finds the final expressions for Qk.LLO, Qk,NLLO and Qk.NNLLO to be: 
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1 + Al''{a)Rks + [Al^'ia)) {Rk.2 - RU) + Al^'{a)Ri, 



These final expression are generic and can be used di- 
rectly. We thus have a procedure which systematically 
sets the PMC scales to all-orders. 

It is easy to see that the leading order values of the ef- 
fective scales are independent of the initial renormaliza- 
tion scale /io- This follows since taking fj-o ^ Q we must 
replace Rk,i — >■ i?fe,i-l-ln(3^/^Q and thus the leading order 
eff'ective scales read In Q^ lo/^o = -Rfc.i+lnQ^/z^O' where 
/io cancels and Eq. (36a) at LLO is recovered. This gen- 
eralizes to any order. Since the /3-function is not known 
to all orders, a higher order residual renormalization- 
scale dependence will enter through the running cou- 
pling. This residual renormalization-scale dependence is 
strongly suppressed in the perturbative regime of the cou- 
pling [37, 38]. 

The effective scales contain all the information of the 
non-conformal parts of the initial pQCD expression for 
p in Eq. (21); this is exactly the purpose of the running 
coupling. The quotient form of Eq. (36c) sums up an 
infinite set of terms related to the known rj^k^Q which 
appear at every higher order due to the special degen- 
eracy of Eq. (21). The method systematically sums up 
all known non-conformal terms, in principle to all-orders, 
but is in practice truncated due to the limited knowledge 



I 

of the /3-function. 

In earlier PMC scale setting [9-11], and its predecessor, 
the Brodsky-Lepage-Mackenzie (BLM) method [7, 8, 14], 
the PMC/BLM scales have been set by using a pertur- 
bative expansion in a and only approximate conformal 
series have been obtained. Here, we have been able to 
obtain the conformal series as revealed in dimensional 
regularization schemes. The final scales in Eq. (36) have 
naturally become functions of the coupling through the 
/3-function, in principle, to all orders. 



B. Automation 

In many cases the coefficients in a pQCD expression 
for an observable are computed numerically, and the {/3i} 
dependence is not known explicitly. It is, however, easy 
to extract the dependence on the number of quark fla- 
vors Nf, since Nf enters analytically in any loop diagram 
computation. To use the systematic method presented in 
this letter one puts the pQCD expression into the form 
of Eq. (21). Due to the special degeneracy in the coeffi- 
cient of the {/3i}-terms, the Nf series can be matched to 



the rj^k coefBcicnts in a unique way^. This ahows one to 
automate the scale setting process algorithmicaUy. 

The n-th order coefficient in pQCD has an expansion 
in Nf which reads: 



Cn,Q + Cn,lNf H h C„^n-lNJ' 



(37) 



By inspection of Eq. (21) it is seen that there are exactly 
as many unknown coefficients in the {/3i}-expansion at 
the order a" as the Nf coefficients, c„j-. This is realized 
due to the special degeneracy found in (21). The rij 
coefficients in Eq. (21) can thus be expressed in terms 
of the c„ J coefficients. The highest power in Nf at any 
order should always be associated with the same power 
in f3o. The ffist /3o appears at order a^~^^. We derive the 
relations between c„j and ri_j for a general gauge group, 
where we define Ca and Cp as the quadratic Casimir 
coefficients of the adjoint and quark representations and 
T as the generator trace normalization. For QCD these 
coefficients read: Ca = N^, Cp = (N^ - l)/2Na and 
T = 1/2. Using that /3o = U/SCa - 4:/3TNf we can find 
r2,o and r2 i: 



r2 = r2fi + npQr2. 



r2.o + r2,i 



UuCa 



An 
r2,i^TNf , 



This leads to: 



3 C2,l 
'"2,1 = -777;^ , r2,0=C2fi 

AI n 



At the next order we have: 



UCa 
AT 



-C2,l 



(38) 



rs = rsfl + n/3ir2,i + {n + l)/^o''3,i + 



n{n + 1) 



/3o^3,: 



where r2,i is already known. Expanding as before in 
terms of Nf (with the higher order /3i coefficient given in 
[29]) we find the matching: 



r3,2 



9 C3,: 



8T^n{n + l) ' 
^3,1 = „,__ ,\,^, [6rc2,i (bCA + 3Cf) 



(n+l)T2 



'"3.0 = C3fi + 



1 



16T2 



-3303,2^^ - 6Tc3,i] , 
[llC^(llc3,2C^+4rc3a) 



(39a) 



(39b) 



-12Tc2,iCa(7Ca + 11Cf)] . (39c) 
Similarly, we can find the r^j coefficients: 



^4,3 



R.2 



rAi 



3! 



AT J n{n + l){n + 2) 



C4,3 , 



32(n+l)(n + 2)T3 



2nT^C2,i (79Ca + 66Cf) - 9 



64(n + 2)r3 



4(3 + 2n) 
n+ 1 



Tc3,2 (5Ca + 3Cf) - 33c4,3Ca - 4Tc4,2 



4r2c2,i (-(37n + 360)CaCf + 2(91n - 150)Cl - 18(n + 6)C|) + 48T2c3,i {5Ca + 3Cf) 



12rc3,2 



+ —CA((l52n + 173)Ca + 33(4n + 5)Cf) - 33C^ (33c4 3CA + STa 2) - A8T^C4, 1 

n+ 1 



ri.o = C4,o + 



64T3 



(40a) 



(40b) 



(40c) 



2T^C2.iCa (8(228 - 77n)C^CF + (840 - 1127n)Cl + 132(n + 6)C|) - A8T'^C3,iCa {7Ca + llCp) 



- 290ATc3,2CaCf + UGT^a^iCA - 1848rc3,2Cl + 484Tc4,2C^ + 1331c4,3Ci 
Using these relations automatically gives the effective scales in Eq. (36c). 



(40d) 



^ In principle, one must treat the Nj terms unrelated to renormal- 
ization of the gauge coupling as part of the conformal coefficient; 
e.g., the Nf terms coming from light- by-light scattering in QED 



and the Nf terms unrelated to the renormalization of the tri- 
gluon and quartic-four-giuon vertices belongs to the conformal 



The automation process can be outlined as follows: 

1. Choose any (5-7?,enornialization scheme and scale. 

2. Compute the physical observable in pQCD and ex- 
tract the A^f coefficients, Ckj- 

3. Find the Pi coefficients, r^^ from the c^j coeffi- 
cients and compute the PMC scales, Qk- 

4. The final pQCD expression for the observable 
reads 



Pfi 



na\{Q) = }Jk+lfia{Qk+l, 



n-\-k 



(41) 



k=0 



This procedure demonstrates that the Nf terms can be 
unambiguously associated to the {/3i}-terms to all orders. 
It also shows that PMC is the underlying principle of 
BLM scale setting. 

The PMC method can be used to set separate scales 
for different skeleton diagrams; this is particularly im- 
portant for multi-scale processes. In general the {ft}" 
coefficients multiply terms involving logarithms in each of 
the invariants [11]. For instance, in the case of qq — > QQ 
near the heavy quark threshold in pQCD, the PMC as- 
signs different scales to the annihilation process and the 
rescattering corrections involving the heavy quarks' rela- 
tive velocity [39]. It also can be used to set the scale for 
the "lensing" gluon-exchange corrections that appear in 
the Sivers, Collins, and Boer-Mulders effects. Moreover, 
for the cases when the process involves several energy re- 
gions; e.g. hard, soft, etc., one may adopt methods such 
as the non-relativistic QCD effective theory (NRQCD) 
[40] and the soft-collinear effective theory (SCET) [41, 42] 
to set the PMC scales; i.e., one first sets the PMC scales 
for the higher energy region, then integrate it out to form 
a lower energy effective theory and sets the PMC scales 
for this softer energy region, etc. In this way one obtains 
different effective PMC scales for each energy region, at 
which all the PMC properties also apply. 

As an important remark, one should keep in mind that 
the determination of the factorization scale is a separate 
issue from renormalization scale setting since it is present 
even in a conformal theory when /3 = 0. Nevertheless, in 
the literature the factorization scale in hadronic processes 
is often set to be equal to the renormalization scale. In 
principle, the factorization scale can be determined if one 
has knowledge of the nonperturbative light-front wave- 
functions of the initial or final state hadrons. However, 
the PMC can also be used to set the scale of the cou- 
pling that appears in the DGLAP or ERBL evolution 
equations, and it is consistent with the usual factoriza- 
tion properties for hard-process cross sections in QCD. 
It is therefore important to separate the renormalization 
and factorization scales in hadronic processes [2]. 



IV. EXAMPLES 



We now consider three examples based on the Adler 
function [43], D, which can be measured indirectly 
through the dispersion relation: 



DiQ') = Q' 



Re 



(s) 



(^ 



ds 



(42) 



where i?e+e- is the ratio for electron-positron annihila- 
tion into hadrons. 

The Adler function is particularly instructive to con- 
sider, since its conformal and non-conformal parts can 
be separated by using RG arguments. Explicitly, the 
Adler function can be written in terms of the vector field 
anomalous dimension, 7, and the vacuum polarization 
function, 11, as follows [44, 45] 



D{Q^) = K-^D{Q^) = 7(a) - /3(a)— n(Q^ a) . (43) 



da 



where /3(a) is the /3- function of the running coupling and 
we have defined the normalized Adler function D where 
K = dp^f Q? and dp is the dimension of the quark 
color representation, which in QCD reads dp ~ Nc- We 
will work with this normalization throughout the related 
examples. In perturbation theory we define 



7(a) = K^7„a(Q)" 

00 

n(a) = K^n„a(0)" 



(44) 
(45) 



which arc now known to four- loop order [45-54]. 
The PMC procedure then follows by absorbing all /3- 
dependent terms, which following Sec. Ill A becomes a 
trivial exercise once the degenerate coefficients r^j have 
been identified. 

As a fourth example, we consider a case where the ex- 
plicit conformal and non-conformal parts are not known. 
Here we make explicit use of the automation procedure to 
derive the special degeneracy as described in Sec. IIIB. 



A. e e — > hadrons 

The ratio for electron-positron annihilation into 
hadrons, Re+e- can inversely to Eq. (42) be computed 
from the Adler function, D, as follows: 



i?e+e-(s) = — 



-s+ze 



D{Q^) 



dQ^ 



(46) 
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It is easy to show that to order a^ the perturbative ex- 
pression for -Rg+g- in terms of 7„ and n„ reads: 

Re+e^ {Q) =70 + iMQ) + [72 + Po^iHQf (47) 

]a{Qf 



2 

+ [73 + AHi + 2/3on2 - /J^"" ^^^^'r^^^ 



[74 + /?2ni + 2/3in2 + 3/?on3 

5 



2-/30/31^ - 3/3o^^ - /3„V^ni]a(Q)^ 



As expected, this expression has exactly the form of 
Eq. (21). with the foHowing identification of the coeffi- 
cients r, oi 



n,t) = 7i 


(48a) 


rj^i ^ Ili-i , i > 2 


(48b) 


?-»,2 = — ^7»-2 , J > 3 


(48c) 


r,,3 = -7r2nj_3 , i > 4 


(48d) 



The expressions for the coefficients 7^ and 11^ can be 
found in Ref. [45, 54], and the four-loops /3-function is 
given in Ref. [29]. The 7^ contain A^y-terms, but since 
they are independent of 5 to any order, they are kept 
fixed in the scale-setting procedure. Notice that this is a 
feature in dimensional regularization. 

Now it is easy to set the exact PMC scales from 
Eq. (36c) using that Rk.j = (-lyrk+jj/jk , 



1.S 



1 ^2 



Ik 

73 

TI2 



72 



1 

2 


da ^ a 




1 
2 


dp , 13 
da ^ a 





-72 



n. 



In- 



Q2 



(49a) 
(49b) 
(49c) 



- (ni + 2 0^ — 71- 31 



/30 + (f)^ 



TT^Hi 



71 -^fn, 



2 \ da 



P da 



71 4\daJ 3 '1 



The final resummed expression for Re+e- reads: 

Re+e- (Q) =70 + 7ia('9i) + 72a(Q2)^ 

+ 730(03)'+ 740(^4)* . (50) 

The scale Q4 is unknown since it requires the knowledge 
of the order a^ coefficient of 11; to leading order it reads: 



In 



Q' 



4,LLO 



Q2 



Hi 

74 



(51) 



however, it is possible to estimate this value. This is so, 
since 114 can be written as 



n4 = - J/?3n2,3 + J/3o'n3,2 + ImoT12,2 

- ^/3on4.i - i/32n2,i - ^/3in3,i + ^, (52) 



where Hij are the coefficients of the bare vacuum polar- 
ization function Ho: 



no(Q,ao) =X!"o 



i-i ( }£_ 

1=1 \^ / k=-l 



J2 f'^^Lk , (53) 



and where e is the dimensional regularization parameter, 
d = 4 — 2e. In Eq. (52) only IIs^o is unknown . We can 
thus compute Qa^llo as a function of IIs^o and for five 
active fiavors we find: 



Q4.LL0 = 0.9 e00ooi3xn.,o g 



(54) 



Because of the small partner it is reasonable to set 
Q4 = Q. The final result in numerical form in terms of 
a = cts/TT for QCD with five active flavors reads: 

Re+e-iQ)^^Re+e-{Q)= (55) 

1 + a(Oi) + 1.84a(Q2)^ - 1.00a(Q3)' - 11.31a(Q4)^ ■ 

This is a more convergent result compared to previous 
estimates, and it is free of any scheme and scale ambigu- 
ities (up to strongly suppressed residual ones). 

From this expression we can determine the asymptotic 
scale A from the empirical data [56]: 

TT-^e+e- (V« = 31-6 GcV) = 1.0527 ± 0.0050 . 

To this end we use the logarithmic expansion for a in 
Eq. (15) and the known expressions for the ji and Hi 
coefficients. For five active flavors we flnd: 



Aa = A 



MS 



419t?68 MeV 



(56) 



which gives us the numerical values for the PMC scales: 
Qi = 1.3 Q ,Q2 — 1-2 Q; Qs — 5.3 Q. These final scales 
determine the effective number of quarks flavors at each 
order of perturbation theory [34] . 

Finally, the QCD coupling at the Mz scale, as{Mz) 
can be computed using again the power expansion for a 
in terms of l/ln(/i/A). We flnd: 



a.(A/z)= 0.132™? 



(57) 



The error on this result is a reflection of the experimen- 
tal uncertainty on R°^_^ which cannot be eliminated. 
This value is somewhat larger than the present world 
average a^iMz) = 0.1184 ± 0.0007, which is a global 
fit of all types of experiments. However, it is consis- 
tent with the values obtained from e~^e~ colliders, i.e. 
as{Mz) ^ 0.13±0.005±0.03 by the CLEO Collaboration 
[57] and as{Mz) = 0.1224 ± 0.0039 from the jet shape 
analysis [58]. Moreover, in computing as{Mz) we have 
assumed masslcss quarks. The estimate will decrease 
when taking threshold effects properly into account"^ as 
shown in [35]. 



We thank Ali N. Khorramian for comments on this point. 
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hadrons 



C. Bjorken and GLS Sum Rules 



It is straightforward to apply our results to the r-dccay 
into hadrons ratio Rr = o-^^i.^+hadrons/o'r^^^+p^+e-, 
which ean be computed from R^+e- [59]: 



RriMr)^^ 



Mt 



Af2 



^^ 1-^ 1 + ^ 



A/2; 



Af2 



Re+e-{s) , 

(58) 



where R^+e- is equal to Re+e~ but with k, = dp^ Q% 
replaced by k! — dF^\Vf f]"^ ., where Vff are the 
Cabbibo-Kobayashi-Maskawa (CKM) matrix elements 
and (^ QfY = 0' since light-by-light scattering does not 
contribute. We define in the same way 7 and 11 (i.e. with 
no light-by-light contributions). In terms of Eq. (21), the 
coefficients for the normalized Rt = Rr/n' read: 



rifl == 7j 



19. 



/265 ^2x i9n^_2 
ri,2 = 1 ^;;; tt I1-2 H , « > 3 



72 



6 



(59a) 
(59b) 

(59c) 



n,3 



265 n\ -- /3355 1977^ 

TT^ ni_3 + 

24 / V 288 12 



7i-3 , J > 4 
(59d) 



The final expression reads 



Rr{Mr) ==70 + lia{Qi) + l2a{Q2f 



(60) 



Since there are three active quark flavors for 
Mr w 1.777 GeV), we find from the CKM matrix 
that K = 3(|KdP + |K.sP) ~ 3. The effective scales read 
Oi = 0.67 Q , Q2 = 0.71 Q, Q3 = 582 Q and for the 
same reason as in the case of Re+e- we set Q4 = Q. The 
scale Qs has been computed to NLLO since its LLO 
value is smaller than the asymptotic scale A. The final 
result in numerical form for three active quark flavors 
read: 



-RriMr) =1 + Q'(Ol) + 2.15a(Q2)' 

+ 3.44a(Q3)^ + 6.64a(Q4)^ 



(61) 



with a == as/iT. Using the asymptotic scale found from 
Re+e~ we estimate the QCD contribution to the r-decay 
to be: 



-1-0.15 



Rr{Mr) =3.66+0:22 



(62) 



This prediction is in good agreement with the ex- 
perimental result from the OPAL collaboration [60]; 
Rf'PiMr) =3.593 ±0.008. 



The Bjorken sum rule and the Gross-Llewellyn Smith 
(GLS) sum rule obey well-known identities in conformal 
field theory, known as the Crewther relations [53, 54, 61- 
64], which through the Adler function can be used to 
expose the conformal terms. In this example, we show 
that both sum rules after PMG scale setting have pertur- 
bative expansions that match exactly the inverse of the 
anomalous dimension, 7"^, and is what one expects in a 
conformal field theory. 

The Bjorken sum rule expresses the integral over the 
spin distributions of quarks inside of the nucleon in terms 
of its axial charge times a coefficient function C^^^: 



rr"(Q2) = f\gr{x,Q')-grix,Q')]dx 
Jo 



9A 

6 



c^^'^(«) + E' 



(63) 



4=2 



where g1^ and gf" are the spin-dependent proton and 
neutron structure functions, gA is the nucleon axial 
charge as measured in neutron /3-decay. The sum in the 
second line of Eq. (63) describes for the nonperturbative 
power corrections (higher twist) which are inaccessible 
for pQCD. Focusing on the perturbative part, we define 

C^J'p(Q2) = 1-3Cf a{Q^) + ^ ^ a{Q^r ■ (64) 

The Gross-Llewellyn Smith (GLS) sum rule. 



GLS/ 



F^P^''''ix,Q')dx = 3NfC'-''^{a), (65) 



relates the lowest moment of the isospin singlet struc- 
ture function F^P+''''{x,Q^) to a coefficient C"^^^(a^), 
which appears in the operator product expansion of the 
axial and vector non-singlet currents. We are again only 
considering the perturbative contribution and define: 

00 
C^i^(g2)=.l-3CFa(g2)-fE C„a(g2)". (66) 

ra=2 

The (extended) Grewther relation [61-64] states that 
there exist a relation between the two sum rules through 
the Adler function D{Q^) given in Eq. (43) as follows: 



5(g2)C^^P(a) = l + ^A» , 
a 

k{a) = aki+ a^ A'2 + a^ A'a -|- 



(67) 



and 



p{a) 



a 
K{a) = a Ki + a^ K2 + a^ K3 



(68) 
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The tilde on D and K indicates the corresponding ex- 
pressions without the Hght-by-Ught type terms, i.e. 



D = D + Diu , 
K = k + Km . 



(69a) 
(69b) 



The term proportional to the /3-function describes the 
deviation from the limit of exact conformal invariance, 
with the deviations starting at order a^. Both sum rules 
have been explicitly computed to four loops and shown 
to obey the extended Crcwther relations [53, 54]*. 

We can use the Crcwther relations to extract the con- 
formal and non-conformal parts of C^^^ and C*^^^. De- 
noting the power expansion of D by 



5(Q2)^l + ^d„a(Q2)", 



n=l 



and expanding its inverse perturbatively gives us 



(70) 



(71) 



C^^s^^-^ ^l^dia + a^ [d\ -d2- (3oKi] 
+ a^ [2did2 -dl^d3 + (3o {diKi - K2) - PiKi] 
+ a* [d\ + dl-di- 3djd2 + 2did3 + /3i (diXi - K2) 
+/3o {-dlKi + diK2 + d2Ki - K3) - I32K1] . 

The expression for C^^^ is the same after putting tildes 
on the coefficients. The di are given in terms of 7^, 11^ 
and /3i as follows: 



di = 71 = 3Cf 



di>2 ^'yi + ^{i-l- k)l3k'ni-i-k 



(72a) 
(72b) 



fe=0 



We use this to find the degenerate r^ j coefficients of 
Eq. (21). 



The degeneracy allows us to resum the series as described 
earlier. The final result is: 

C«^^(a) ^ 1 - a(Qi)7i + a{Q2f (7? - 72) 

+ aiQsf {-jf + 27271 - 73) 

+ a{Q4)^ {jt - 3727? + 27371 + 7I - 74) + O (a') , 

(74) 

exposing the r^ coefficients. This expression is simply 
the inverse of the anomalous dimension: 



C^^^(a)==7"'(Qi,Q2,g3,...) 



(75) 



where we used that 70 = 1. The arguments of 7"^ 
on the right hand side indicate the effective scales 
at each order in perturbation theory, once the in- 
verse is Taylor expanded. All the above expressions 
also apply to the Bjorken sum rules, with the coeffi- 
cients replaced by the ones with tilde. In particular, 

C^^P{a) = rHQuQ2,Q3,---)- 

Since, the Adler function itself after PMC scale setting 
is simply given by the anomalous dimension: 



D{Q)=i{Qi,Q2,Q3,- 



(76) 



and correspondingly for D, the Crcwther relations can 
be expressed as 



lit^i, ^^2, ■ ■ ■) 

7(M1,M2,---) 



(77a) 



(77b) 



where the last equality follows due to conformality. These 
are the generalized Crcwther relations, which set the 
commensurate scale relations between the scale of the 
Adler function and those of the sum rules. 



r2,i = -K^ - Hi , 
r3,i--^-n2+f^+ni)7i 



''4,1 



2 
K3 



71 



r4,2 = -(^ini + n2), 

'"3,2 = , r4.3 = . 



(73a) 
(73b) 



--n3 + (/i2 + 4n2, g 

-(':|i+ni')7f + (Xi + 2ni)|, (73c) 



(73d) 
(73e) 



* There is a recent claim [55] that the existing four-loop coefficient 
of the Bjorken sum rule [53, 54] is missing some singlot-diagram 
contributions. This is relevant only for the explicit evaluation of 
K3, and does not change the results of this section. 



D. Static Quark Potential 

As a last example we consider the potential between 
two static quarks, where the degeneracy is not explicitly 
apparent in the literature. The static quark potential is 
known to order a* in the MS-scheme as an expansion in 
the number of massless flavors, Nf [65-70]: 



ViQ'') 



(4^)^Cj 



a{Q'') 1 + (c2,o + C2,iAf/)a(Q') 



+ (c3,o + C3,iA^/ H- C3,2Nf)a{Qy 

+ (C4,0 + C4,lA^/ + C4,2iV/ + C4,3Af/)a((3^)^ 



+8^2ciln^a(Q2)3 



+ Oia'), 



(78) 



where we have chosen the initial scale of the running 
coupling /i^ = Q^, but kept the explicit IR divergent 
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find the degenerate coefficients of tlie static potential to 
be: 



logaritlim Infi'jj^/Q'^, wliich is not related to coupling the explicit expression for Cij given in Refs. [68-70] we 
constant renormalization, but is coming from the non- 
abelian gluon 'H-diagram'. This IR divergence is a fea- 
ture of pQCD and is canceled by non-perturbative contri- 
butions from the 'ultra-soft' region [70, 71], which is con- 
trolled by the domain of color-confinement. The regular- 
ization comes from the energy difference between color- 
singlet and octet intermediate states [71]. 



'2,0 



:Ca 



o V _ 

"q ' ''3,2 — 



' '4,3 



The degenerate coefficients rtj are determined from 
the Cij-coefficients as given in Sec. IIIB for n ~ 1. From 



n-\-l,n 



(79a) 



(79b) 



'^^0 =^ [(532 - 1584C3 - 97r4 + lUn^) Ca + 33 (48(3 - 35) Cp] 



'3,1 



7C3-f la 



35 



- 6C3 Cf , 



11 / ^abcd yjabcd 

rlo =7n:CA ( (456C3 - 1440C5 + 571) C| - 9 • 56.83(1)^ ^ 



36 



TNa 



/3077C3 



4847r^ 



Ca -^7^ - 1293.54(1) + — — - 136.39(12) ^ ^ 
V 3 135 



jabcd ^abcd 



V _/^ 392(3 



'4,2 



4OC5 
23(3 + 



66769\ 
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38C3 571 



r; , = I — - 2OC5 - ^^1 CaCf + C\ ( 196.58 - 192(3 - 



Na 



135 



36 



)c| + 56.83(1)^ 



'abcdjabcd 



47r4 2981 \ ^ /5171 



ATNa 
- 26C3) Cf . 



758C3 
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The static potential can then be written without any explicit dependence on Nf: 



V{Q')=- 



a{Q') 1 + (<o + rl^MaiQ') + (rg^ + (i^rl^ + 2f3^Tl^ + litrlaXQ') 



3709\ 2^ 



+ «o + p2rX,i + 2P^rl, + ^/3i/3or|;2 + 3/30^1 + 3/?o'<2 + P^MQ^ +&ti^C\ In ^a{Q') 



2\3 



(80a) 
(80b) 

(80c) 

(80d) 
(80e) 



(81) 
+ 0(a5) . 



r 



Next, we reduce the expression to the conformal series by 
using the PMC scales, which are read of from Eq. (36): 



V{Q^) = - 



{A^YCf 



aiQD+rl^aiQlf+rlMQlf 



2\i 



-rAoa{Qi) 



Sn'Clln^aiQr 



(82) 



gauge theory with Nf light quarks read: 

ay(Q2) ^a(Qi) - 0.64a(Q2)' - 0.78a(Q3)^ 



3.49 



27 ^,'' 
8^^"q^ 



c^iQir 



with 

Ql 



all orders 
Q2,NLLO 



.Q2exp(-5/3) , 

=Q2 exp (0.42 4- 0.57/?oa((3^)) 



3.LLO=Q'exp(-5.87), 



(83) 

(84a) 
(84b) 
(84c) 



The expression in the bracket defines the effective charge 
ay^Q^). Its IR divergence can be removed by adding the 
ultra-soft contributions to the static potential. The final 
expression for the effective charge ay = Anay in a 5*1/(3) 



where ^0 = H - f/V/. Note that the PMC scale for Qi 
holds to all orders, which follows since r„+i_„ = (5/3)" 
at any order n. The scales for Q2 and Q^ have been ex- 
panded in a to the order consistent with the pQCD trun- 
cation. We note that the third PMC scale Q3 is greatly 
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suppressed compared to the kinematic scale Q. This is to 
be expected physically, since the kinematically accessible 
region shrinks with the loop order. This indicates the 
breakdown of perturbation theory in higher order QCD, 
which is also explicitly evident by the IR divergent term 
appearing in the a^ coefficient. At these higher orders 
non-perturbative effects must be taken into account. 

Finally, we can now explicitly show that PMC scale 
setting is consistent with the Cell Mann-Low (GM-L) 
scheme and the effective coupling in QED. It is well 
known that the effective QED coupling in the massless 
limit, is related to the MS coupling by a scale displace- 
ment; aGM-L(Q^) = oiyf^iQ'^e^^/^). The effective QED 
coupling is precisely deiined as the effective charge of the 
QED static potential between two (formally) infinitely 
charged particles. By inspection of Eq. (79a), (80a) and 
(80c) it can be seen all the higher order conformal co- 
efficients, ri>2,o are proportional to non-Abelian group 
invariants, which vanish in the Abelian limit [36]; e.g. 
Ca -> 0. This generalizes to any order. This means that 
the PMC expression for the effective charge of the static 
potential is given to all orders in perturbation theory by: 



av,QED{Q^) = a{Q^ 



-5/3 



(85) 



Thus PMC reproduces the correct result in the Abelian 
limit. 



V. COMMENSURATE SCALE RELATIONS 



where aiz = ria-jz, the effective charge itself can be con- 
sidered as a running coupling of a physical scheme re- 
lated to the corresponding observable. The above expan- 
sion then defines the scheme transformation from the TZs- 
scheme to the p-scheme. Since any two effective charges 
ga and as can be computed in the TZs-scheme, it fol- 
lows that aA can be written as an expansion in as by 
scheme transformations. Thus, any effective charge de- 
fines a physical renormalization scheme. The /3-functions 
of such schemes are different from the /3-function of the 
7?.a-schemes, but are related by the identity: 

M^a) = ^M&n) , (88) 

where A is some physical scheme corresponding to the 

effective charge a a- From this identity it follows that the 

first two coefficients of the /3-function are universal [72] . 

The expansion of ua in a-jz can be put to the form: 

aA{Q'') = aniQ^) + [4,0 + ^ori,]aniQ^f 

+ Ho + ^irt + Worti + $lri^2]&niQ'f 
+ Ho + Pfri, + 2p,ri, + ^PiPori^ + iPort 
+ iPlrl^ + PlrUan{Q^f + 0{a\) , (89) 

where rf, are related to the coefficients r^ ., of the ob- 
servable A, by: 



We demonstrate that the generic expression in Eq. (21) 
extends to any scheme, that is, the special degeneracy in 
the TZs-scheme of an observable is inherited in all physical 
schemes. This is done by relating different observables 
in pQCD using the effective charge method [4-7]. These 
commensurate scale relations must be independent of the 
choice of scheme. The scales are given by the systematic 
scale-setting method just described. 

Any observable p can be used to define an effective 
charge Op. Considering the case where the Born level 
result for the observable is just a constant such as Re+e- 1 
i.e. n = in Eq. (16), the effective charge is defined by 
the relation 



p(Q') = Po(Q')[l + ap(Q' 



(86) 



where po is the Born (tree-level) result and Q^ is the mea- 
sured scale. Thus, Up can be understood in perturbation 
theory as summing up the entire perturbative series into 
one effective coupling; the effective charge of the process. 
It follows that the effective charge has an expansion in 
the TZs coupling a-jz similar to the expansion in Eq. (21). 
(we put back the index TZ on the coupling in this section 
to avoid confusion). By normalizing running coupling 
such that at leading order the running coupling is equal 
to the effective charge; i.e. 



ap = a-jz 



r2.2 

72 ■R- 



^3.3 



(87) 



A ^ ''J 
» j' i-j 



(90) 



and $i are here the coefficients of the beta function of 
a-ji, which are related to the usual /Si in Eq. (5) by: 



A. 



(91) 



The renormalization scheme dependence of /32 is denoted 
with a superscript. The coefficient /?2 can be found in 
terms of /3o, Pi and /3^ from Eq. (89) using Eq. 



(3^=^f-^4, + ^^,{ri^2 



'2.1 



+ 2/3o^ {rt - riort) " Po{4,o + ^3^0) ■ (92) 

Using this, it can be shown that the special degeneracy 
in Eq. (89) is preserved when relating the effective charge 
a A with another effective charge as, i-C. 



„AB 



aA{Q') = asiQ') + <S*+/3o<f aeiQ') 



„AB 



AB 



q2AB 



riS+Pirt?+2f^ort?+P'ori 



aBiQ') 



(93) 



2\3 



AB 



r^o + f^"rq- + 2/3ir3-f + -fiifiorta 



■3/3o' 



,AB 

4,1 



m 



2„AB 

4,2 



Pirti asiQY 
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where the coefficients rff are related to the TZs- 
coefficients as follows: 



*j 



'2,0 

'^2,1 

^AB 
'3,0 

^AB 
'^3,1 

^AB 
'3,2 

^AB 
U,0 



„AB 



„AB 



^AB 
4,3 



'2,0 



'2,1 



'3,0 



'3,1 
'3,2 



'4,0 



'2,0 I 
'^2.1 ' 



2r. 



S :„AB 



„B AS 



'3,1 
'3,2 

'^4'!o - 3r. 



'2,l'2,l 
2,0'3,0 



R AB 



'2,0 



(94a) 
(94b) 
(94c) 
(94d) 
(94e) 



2r: 



i3 \^AB 
3,0J'2,0 



(94f) 



O^S AS 
^'2,0'3,1 



2r. 



'2,l'3,0 



B „AB 
2,l'3,0 

.5 ^-B \„AB 



„A B 
/on'' 



'4.2 



r 



4,2 



l/j.-B , 9 « H ■. AB 
ol'3,l + ^'2,0' 2,17' 2,0 1 



3,l'2,l 



^ '"2,0'"3,2 ^ ■o(''2,l + 2^3 2)^2 



'4,3 



'4,3 •J'2,1'3,2 



«"3,2'2,0 



3,0' 2,1 

(94g) 

(94h) 
(94i) 



This demonstrates that the special degeneracy of the {/3i} 
coefficients is not a prerogative of the 7?,5-schenics, but is 
a general feature of perturbation theory. Since any two 
effective charges are related by the same perturbative 
pattern of Eq. (21), we can directly use the systematic 
scale-setting method presented in the previous section 
to eliminate the initial scale ambiguity in Eq. (93). This 
explicitly shows the renormalization-seheme invariance of 
the scale-setting method to all orders in perturbation the- 
ory. The final relation between any two effective charges 
is thus: 



' 3,o''«B(g^)^ + rt^^aBiQir + 0{a%) , 



„AB, 



-)2\3 I „AB , 



(95) 



where the commensurate scale relations between the two 
charges are exactly the PMC scales Qi as given by 
Eq. (36) with the /3-function being that of ub ■ For com- 
pleteness we provide also the expression for the four-loop 
/3-funetion coefficient of an effective charge A: 



'3,2 



+ 4/3i/3o (r^i - '^^o'^^i) + /3i<o 
+ 6/33 f 2r^nr^r - 2?- 



'4,3 



''3,1^2.1 



^A J\ 
'2,0'3,2 



^A 

'4.2 



+ 



6/3o' (2r2^i 
2/3o (2r, 



2,0 



2ri 



A 3 
2,0 



3''3,0''2.0 



'4,0 



(96) 



As a particularly simple example, we relate the effec- 
tive charge of Rt] a,-, to that of Re+e- ! O'R a-^d apply the 
systematic scale-setting method to derive commensurate 
scale relations between the two effective charges. The 



final result is completely independent of the intermedi- 
ate renormalization scheme and scale used to compute a,- 
and ur. 

The degenerate coefficients of Or and afl in the TZs- 
scheme can be read of from Eq. (59) and (48), from which 
we compute rp- . Using Eq. (94) and Eq. (95) we can 
readily express Ur as a perturbative series in gr: 



ariQ^) ==ai?((9fl,i) 



r{ 



-MQI..3)' 



7174. /6Z — 37273, ;fc; 



MQIa)^ , (97) 



where the PMC scales Qn^i are given by the systematic 
method and where "fi^wi is the light-by- light part of 7^; 
i.e. the two effective charges are equivalent up to light- 
by-light terms. The corrections start at order a^, since 



there is no light- by-light contribution to 72. 
scales expanded in ar{Q) read: 



The PMC 



In 



Q' 



R,l 



In 



Qrs 



19 169 
12 " 144 

f 761/3g 
l^ 192 

19 

12 



169^1 \ 
96 ] 



2\2 



aniQ') 



ni \ 73,;fc; ^3dbi 

71 



73 



73 



(98a) 
(98b) 



By definition the scale where the expression for the 
effective charge Qt applies is Q^ = M^. At this 
scale, the number of light flavors is Nf = 3. Light- 
by-light diagrams are proportional to (^fQf)^, which 
vanishes exactly when summing over the three light 
quarks. Therefore in three-flavor QCD the two ef- 
fective charges are identical to all orders; i.e. using 
flfl/r = 7iai?/r/(47r) = aR/r/Tr, we have 



a,(M2) _ «fl(Q|j,i) 



(99) 



where the commensurate scale for ur up to four-loop 
order is given by: 



In 



^R,l 



19 169 ai?(M2) 83273 afl(Af 2)2 
12 



64 



3072 



(100) 



This relation (at one lower order) has been shown to be in 
very good agreement with experiment [73] demonstrating 
a highly nontrivial consistency check of QCD, free of any 
scheme and scale ambiguities. 



VI. CONCLUSION 

In this paper we have shown that a generalization of 
the conventional MS-schcme is illuminating. It enables 
one to determine the general (and degenerate) pattern 



16 



of nonconformal {/3i}-terms and to systematically deter- 
mine the argument of the running coupling order by or- 
der in pQCD, in a way which is readily automatized. 
The resummed series matches the conformal series, in 
which no factorially divergent n!/3"a" "renormalon" se- 
ries appear and which is free of any scheme and scale 
ambiguities. Thus using the PMC/BLM procedure, all 
non-conformal contributions in the perturbative expan- 
sion series are summed into the running coupling by shift- 
ing the renormalization scale in a^ from its initial value, 
and one obtains unique, scale-fixed, scheme-independent 
predictions at any finite order. The resulting PMC scales 
and finite-order PMC predictions are both to high ac- 
curacy independent of the choice of initial renormaliza- 
tion scale. The PMC procedure also provides scale-fixed, 
scheme-independent commensurate scale relations, rela- 
tions between observables which are based on the un- 
derlying conformal behavior of QCD such as the Gener- 
alized Crewther relation. Furthermore, we have shown 
that PMC is consistent with QED scale-setting, where 



there is no ambiguity in choosing the final scale of the ef- 
fective coupling. The PMC satisfies all of the principles 
of the renormalization group: reflectivity, symmetry, and 
transitivity, and it thus eliminates an unnecessary source 
of systematic error in pQCD predictions. 
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